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ABSTRACT 



We report thorium abundances for 77 metal-poor stars in the metallicity range of —3.5 < [Fe/H] < —1.0, based on "snapshot" 
spectra obtained with VLT-UT2/UVES during the HERES Survey. We were able to determine the thorium abundances with better 
than 1(7 confidence for 17 stars, while for 60 stars we derived upper limits. For five stars common with previous studies, our results 
were in good agreement with the literature results. The thorium abundances span a wide range of about 4.0 dex, and scatter exists in 
the distribution of log(Th/Eu) ratios for lower metallicity stars, supporting previous studies suggesting the r-process is not universal. 
We derived ages from the log(Th/Eu) ratios for 12 stars, resulting in large scattered ages, and two stars with significant enhancement 
of Th relative to Eu are found, indicating the "actinide boost" does not seem to be a rare phenomenon and thus highlighting the risk 
in using log(Th/Eu) to derive stellar ages. 
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1. Introduction 

The observed chemical abundances of metal-poor stars provide 
rich information on star formation and nucleosynthesis in the 
early stages of evolution of our galaxy and other galaxies. In 
particular, the study of heavy elements beyond the iron group 
in metal-poor stars have greatly improved the understanding 
of the neutron-capture nucleosynthesis processes in the early 
Universe, e.g., placing constraints on the astrophysical site(s) 
of r-process. Neutron-capture elements in r-process enhanced 
metal-poor stars have been reported to match the scaled solar 
r-process abundance pattern a t least in the range of Z — 52 ^ 76 
(see e.g., ICowan et alj l2002t ISneden et alj 120031: iHonda et alj 
12004 *). In contrast, lighter elements (Z < 50) abundances devi- 
ate significantly from the Solar-pattern, indicating the existence 
of different in itial producti on ratios for the r-process elements 
dSneden et aLi2000, 2003: Honda et al.l2004 . In old, metal-poor 
stars, thorium and uranium with relatively long half-life are the 
only heavier r-process elements we could observe today, and of- 
fer a way to bette r understand the heav iest products in r-process 
nucleosynthesis dRoederer et al.ll2009l) . 

The detection of the radioactive-decay elements in r-process 
enhanced metal-poor stars allowed a new approach to derive the 
age of the oldest stars, by means of comparing the observed 
ratios of radioactive decay elements over stable elements with 
the corresponding initial values at the time when the star was 



borrQ, hence providing a lower limit on the age of the uni- 
verse. The actinide thorium is one such radioactive species. It 
is produced in the r-process, and has been widely used, par- 
ticularly log(Th/Eu), for cosmochronometry to es timate stel- 
lar ages (e .g., Snedenetal. 1996; Hill et al. 2002; Honda et al.l 
2004 lF^bel et al.ll2007l: IHavek et al.l ii2009). Some authors have 
argued that log(Th/Eu) is not a reliable chronometer to de- 
rive stell ar ages, because of the r elatively large uncertain- 
ties (e.g.. iG orielv & Clerbaux"1999^, lOorielv & Arnouldl 120011; 
iHonda et al.l Eb04; Frebel et al. 2007). Due to its longer half-life 
of 14.05 Gyr as compared with, say, 4.5 Gyr of uranium, some 
authors have suggested that the log(U/Th) pair would be a better 
chronometer, because of their close atomic n umber and thus rel- 
ativ ely small nuclear physics uncertainties dCavrel et aLlbOOU 
Wan aio et al] [2002; Beers & Christlieb 200'3. Additionally, be- 
cause Pb and Bi mostly originate from a-decay of the Th and 
U isotopes, the measurement of Pb or Bi abundances can of- 
fer a con sistency check on t he calculated initial abundances for 
this pair dCowan et al.lll999l) . However, it is very difficult to de- 
tect U lines in stellar spectra, due to the blending of the weak 
U lines (e.g., 3859 A ) and the relatively low amo unt of ura- 
nium in an old star with ag e of ^ 12 — 15 Gyr (see iPIez et al.l 
2004; R oederer et al.ll2069l) . Alternatively, Hf behaves similarly 
to third-peak elements and the log(Th/Hf) ratio has been sug- 
gested bv Kratz et al. (2007.) as a promising tool for chronometer 
studies. 



* Based on observations collected at the European Southern 
Observatory, Paranal, Chile (Proposal Number 68.B-0320). 



Accurately speaking, it was the time when the radioactive decay 
elements were produced, not necessarily to be the same as when the 
star was born. 
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Many efforts have been made to determine thorium abun- 
dances in metal-poor stars using high resolution spectra. So far, 
more than 30 metal-poor stars ([Fe/H] < —1) have been re- 
ported with thorium abundances in previous studies, and the 
error is typically of order of 0.15 dex. It appears that among 
the r-process enhanced stars with measured Th or other actinide 
elements (e.g.,U), except for an actinide normal group, which 
means no obvious enhancement of actinide element abundances 
with respect to the scaled solar r-process pattern, an actinide 
boost group (e.g., CS31082-001, HE1219-0312, CS30306-132) 
also came to be known, although the former group seems to 
be more common among r-process enhanced metal-poor stars 
(e.g., CS22892-052, CS29497-004). This imples that for ele- 
ments in the range of Z > 90, significantly different chemi- 
cal yields might be produced due to the various conditions of 
the star formation region s ([Hill et al. 2002; Roederer et al. 200^ 
iMashonkina et al.ll2010h . 

To obtain accurate age estimates and to explore r-process 
nucleosynthesis in metal-poor stars in detail, high-precision Th 
abundances are needed. Further, large samples of stars with Th 
abundances are also important to study the distribution of tho- 
rium abundances in metal-poor stars. The HERES survey of- 
fers a good opportunity to perform such a study. More than 22 
elemental abundances (not including Th) for the sample stars 
have been reported in previously published papers in this se- 
ries(e.g., Barklem et al. 2005; Zhang et al. 2010), and thorium 
abundances have been studied furtherly with be tter qual ity spec- 
tra in a few stars tHavek et al.ii2009i: iMashonkina et al.ll2010i) . 

Here we report thorium abu ndances for 77 stars from the 
sample of iBarklem et al.l (|2005') (hereafter Paper II), using a 
modified version of the analysis method described in that paper. 
A brief description of the sample is given in sect. [2] The abun- 
dance analysis is described in sect. [3] In sect. |4l thorium abun- 
dances are presented and compared with previous works, and 
the implications for stellar age estimates are discussed. Section 
|5]presents the conclusions. 

2. The sample 

This work is based on the moderately high-resolution "snap- 
shot" spectra of 253 HERES stars. Readers are referred to 
IChristlieb et al.l (l2004h and Paper II for detailed information 
on the observations and the sample selection. For convenience, 
we repeat the most important basic information on the sample. 
The spectra (R - 20000, A = 3760 - 4980 A , and typical 
S/N ~ 30 to 50) were obtained during the HERES survey with 
ESO-VLT2/UVES. A total of 373 spectra were observed, which 
was reduced to 253 stars in Paper II when stars were removed 
for various reasons, most importantly due to strong molecular 
carbon features leading to significant blending. Much of the 
remainder of the sample has been analysed by iLucatello et al.l 
(i2006). These 253 stars are the starting point for our analysis; 
however, reasonable estimates of Th abundances or even upper 
limits could not be obtained for all stars for reasons that will be 
discussed further below. 



3. Abundance analysis 

We derive abundances from the Thll 4019.12 A line, the only 
Th line strong enough to be detected in our spectra. This line 
is unfortunately blended, including blends of ^'^CH. The line 
has been analys ed using the automated sp ectrum analysis code 
based on SME dValenti & Piskunovll 19961) described in Paper II 



Table 1. Line list. 



Species 


A 




log gf 


Refs. 




[A] 


[eV] 






Th II 


41)19.129 


0.000 


—0.228 


1 


CH I 


4019.000 


0.460 


-1.163 


2 


^^CH I 


4019.170 


0.460 


-1.137 


2 


Pel 


4019.050 


2.608 


-2.780 


3 


Col 


4019.110 


2.280 


-3.287 


2 


Col 


4019.118 


2.280 


-3.173 




Col 


4019.120 


2.280 


-3.876 




Col 


4019.125 


2.280 


-3.298 




Col 


4019.125 


2.280 


-3.492 




Col 


4019.134 


2.280 


-3.287 




Col 


4019.135 


2.280 


-3.474 




Col 


4019.138 


2.280 


-3.173 




Col 


4019.140 


2.280 


-3.298 




Col 


4019.272 


0.580 


-3.480 


2 


Col 


4019.281 


0.580 


-3.470 




Col 


4019.294 


0.580 


-3.220 




Col 


4019.296 


0.580 


-3.330 




Co I 


4019.322 


0.580 


—4.090 




Col 


4019.332 


0.580 


-4.040 




Nil 


4019.058 


1.935 


-3.174 


2 


Ce II 


4019.057 


1.014 


0.093 


2 


Nd II 


4018.836 


0.060 


-0.880 


2 



Ref erences. (1) - iNilsson et all ( l2002h . (2) - iJohnson & B"oit3 j2001h 
(3) - iKupka et all n999l) 

and used in lJonsell et al.l (l2006l) (hereafter Paper III). The stellar 
parameters and abundances of other elements, particularly those 
giving rise to blends with the Th line, were adopted directly from 
Paper II. In order to be able to model the ^''CH blends, some 
modifications were made to the code to enable estimates of the 
^^C/ ^^C ratio to be obtained. This ratio has been determined 
from isolated ^^CH features between 4210 and 4225 A and at 
4370 A, as shown in Paper III, though now in an automated man- 
ner 

Th e employed line data is shown in TablefT] and is essentially 
that oflJ ohnson & Boltd (12001 ). The / value for the Th II line has 
been updated (a change of -1-0.05 dex) and the wavelengths and 
excitation potentials have been ne gligibly changed in some cases 
based on VALD values dKupka et al. 1999). All blends within 
0.5 A either side of the Th line in the list of lJohnson & Bolt3 
(1200 II) are included, except a line of U and and a line of V, which 
we are unable to model since we do not have abundances for 
these elements; their contributions are expected to be negligible 
in any case. A window 0.3 A wide centred on the Th line is used 
for the fitting. 

In Paper II we required at least one line to be detected at 
the 3(7 confidence level to claim an elemental detection. Such 
a high threshold would result in a very small number of detec- 
tions given the weakness of the single Th line we used and the 
quality of our spectra. Here we reduce this requirement to la, 
though with some additional constraints. First, the significance 
of detection was caculated from the fit without blending, i.e., 
m' = m X (1 — db/dfit), where db and dgt are the line depths 
of the blendings (without thorium line) and the whole fit respec- 
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lively; m is the detection level of the the whole fit. Second, all 
lines calculated to have m above the Icr detection level were 
subjected to manual inspection and adjustment adopted for some 
stars, usually due to noise or severe blending casting significant 
doubt on its validity. Finally, except for several stars raised to de- 
tection level considering their strong line strengths, for all other 
stars with detections rn' below the 1 a confidence level and those 
rejected during manual inspection, we calculated Icr upper lim- 
its, but rejecting those with uncertain upper limits, i.e., blend- 
ings occupying more than 80% of the fit or one of the fits not 
converged. 

The errors in the abundances could in principle be estimated 
using the methods of Paper II, but is complicated by the addi- 
tional uncertainties from blending. In Paper II we found rela- 
tive and absolute errors of ^ 0.18 and ^ 0.25 dex, respectively, 
for species with similar atomic structure to Th II, and thus lines 
with similar sensitivities to stellar parameters, e.g. Nd II, Sm II, 
Eu II. These cases have more and stronger lines and no signifi- 
cant blending issues. Accounting for these extra sources of error 
we estimate relative and absolute errors in the Th abundances of 
~ 0.25 and ~ 0.3 dex, respectively. Similarly, by comparison 
with errors in Nd/Eu etc., an error of ^ 0.25 dex is estimated for 
the abundance ratio log(Th/Eu). 

However, it should be noted that, in our analysis, we have 
adopted the detection significance of Icr, which means the statis- 
tic uncertainty will dominate the general error. Thus, we value 
any interesting results from our analysis, but are very cautious to 
give any conclusion based on our results. 



4. Results 

The detected thorium abundances for stars classified into differ- 
ent subclasses are tabulated in Table |4] In Table |5] (online only), 
stars with upper limits are included, together with the atmo- 
spheric parameters and some other abundances of interest from 
Paper II. The confidence level of the detection ma and m'cr, 
computed as described in section 3.4 of Paper II and in section 
[3]of this paper, are also tabulated. The errors are not listed since 
they are either the same for all stars in the cases involving Th 
(see above), or available in Paper II. 

Throughout this discussion w e classify different types of 
neutron-capture stars according to lBeers & Christliebl (l2005h as 
summarized in Table [3] In addition, we adopt the definition of 
s-II stars from Paper II. Compared with the s subclass, s-II stars 
actually are those with both strong r- and s- process enhancement 
([Ba/Fe] > 1.5, [Eu/Fe] > 1.0 and [Ba/Eu] > 0.5). 

We note, HE122 1-1948 is one of the several C-enhanced 
stars ([C/Fe1=l . 42) in our analysis, which has been analysed by 
iLucatello et al.l (l2006l) manually, giving a metallicity of —2.60 
and [Eu/Fe]=2.1. Since in Paper II, it was mentioned that our 
method is not suitable for C-rich stars spectra, and actually this 
star was noticed in figure 5 of Paper II with A[Fc/H] > 1 com- 
paring the final [Fe/H] with the initial estimate, thus in this work 
we demote this star to an upper limit, although m' > 1 was ob- 
tained, which we subjected to an overestimated S/N. HE1221- 
1948 is the only carbon rich one ([C/Fe] > 1) among the most 
metal-poor stars ([Fe/H] < —3.0) with detected thorium abun- 
dances in our sample. 



4. 1. Thorium abundances distribution 

In the sample of 253 stars, finally we obtain thorium abun- 
dancefl for 17 stars. For another 60 stars, we obtained upper 
limits. In Fig. [1] we give the thorium abundances distribution 
of these stars. We were unable to derive the thorium abundance 
for the other 176 stars, 79 of which were set aside for various 
reasons, the most common being that the spectrum was too con- 
taminated with blends to derive anything reliable given the reso- 
lution of the observed spectrum. We adopted Icr detection level, 
and since our aim is to investigate the general trend of thorium 
abundances, this detection level is adequate to give a reliable 
distribution. 

As shown in FiglT] in the distribution of thorium abundances 
by the Icr detection, majority of the stars lie in the range of 
—2.0 to —1.0, about 82% of all the detections, and the distri- 
bution including Icr upper limits extends from —3.0 to 1.0, hav- 
ing a peak at —2.0 to —1.0. The thorium abun dance results in 
previous studies all fall in this range (see e.g., Roeder er et al.l 
2009; Honda et al. 200 4). In the lower panel of FigH] the only 
one detection above 3cr is owned by CS31082-001, which has 
loge(Th) = —1.00, and including those with 3cr upper lim- 
its the distribution spaning from —3.0 to 1.5, has a peak at the 
range of —1.0 to —0.5. Only the higher Th abundances side of 
this histogram can give useful information about the true shape 
of the cosmic distriubtion, since there is a natural detection bias 
towards larger abundances. The existence of a large scatter in 
thorium abundances among metal-poor stars as seen from this 
histogram can be explained by the poor mixing in the early uni- 
verse. 

However, because the abundance of thorium is based on one 
single Th II line 4019 A , it should be borne in mind the that Th 
abundance is susceptible to errors due the difficulties in mod- 
elling the blends. The line Th II 4019 A is blended with tem- 
perature sensitive spectral features, such as Co and Ni, thus, we 
expect some degree of the scatter is due to errors in the analysis, 
e.g., uncertain Teff, rather than real cosmic scatter However, as 
shown in FiglT] and Fig.[3]of this paper, the scatter significantly 
exceeds the estimated errors, and thus given the detection bias, 
it gives a lower limit to the range of the true cosmic distribution. 

4.2. Comparation witli previous studies 

As shown in Table |2] we obtained Th abundances for 
five stars which have been studied previously in detail 
with similar or better quality observational data by others, 
namely CS 22892-052, CS 29497-004, CS 31082-001, 
CS 29491-069 and HE 2327-5642. Note that we failed to 
derive the thorium abundance for a well studied bright halo 
star HE 221170, since our modelling indicated Th II line in 
this star to be severely blended by Co lines; th is star has 
been reported thorium enhance d compared with iron divans et alj 
l2006l:lYushchenko et alj|2005l) . The same for another strongly r- 
process enhanced star HE 1219—0312, in which the blendings 
exceed 80% of the line fit, resulting a very uncertain abundance. 
HE 0338—3945 is an reported s-II carbon enhanced star ([Eu/Fe 
= 1.89, [Ba/Eu] = 0.52, and [C/Fe] = 2.07]), which experienced 
strongly both r- and s- process enhancements. For the same rea- 
son as for other carbon-enhanced stars, our method failed to get 
a reliable thorium abundance for this hot dwarf. For rest of the 



" Here we adopt the standard spectroscopic notations: 
loge(A) = logio(iVA/A'^H) + 12 for abundances; 
[A/B] = logjo(iVA/A''B)*-logiQ(A''A/A^B)0 forrelative abundances; 
log(A/B)=loge(A)-loge(B) for abundance ratios. 
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Table 2. Comparison with literature results for seven common stars 









logg 


[Fc/H] 




loge(Th) 


Notes 






[K] 


r — 2n 

[cm s J 




r. —In 

[kms ^] 






CS 22892 


-052 


4884 


1.81 


-2.95 


1.67 


-1.76 


This work 






4790 


1.6 


-2.92 


1.8 


-1.42 ±0.15 


Honda et al. (2004) 






4800 


1.5 


-3.1 


1.95 


-1.57 ±0.10 


Sneden et al. (2003) 






4710 


1.5 


-3.2 


2.1 


-1.60 ± 0.07 


Sneden et al. (2000) 






4800 


1.5 


-3.1 


1.95 


-1.60 ±0.13 


Roederer et al. (2009) 


CS 29497 


-004 


5013 


2.23 


-2.81 


1.62 


-1.17 


This work 






5090 


2.4 


-2.81 


1.6 


-0.96 ± 0.15 


Christlieb et al. (2004) 


CS 31082 


-001 


4922 


1.90 


-2.78 


1.88 


-1.00 


This work 






4825 


1.5 


-2.9 


1.8 


-0.98 ± 0.05 


Hill et al. (2002) 










-2.9 




-0.98 ±0.05 


Plez et al. (2004) 






4790 


1.8 


-2.81 


1.9 


-0.92 ±0.10 


Honda et al. (2004) 


CS 29491 


-069 


5103 


2.45 


-2.81 


1.54 


< -1.37 


This work 






5300 


2.8 


-2.6 


1.6 


-1.46 ±0.25 


Roederer et al. (2009) 






5300 


2.8 


-2.6 


1.6 


-1.43 ± 0.22 


Havek et al. (2009) 


HE 2327- 


5642 


5048 


2.22 


-2.95 


1.69 


-1.45 


This work 






5050 


2.34 


-2.78 


1.8 


-1.67 ±0.20 


Mashonkina et al. (2010) 



15 



10 





25 

20 

15 
10 




1-a 



-4 -3 



-2 -1 
loge(Th) 



3 -a 



Fig. 1. Histogram of derived thorium abundances. Icr (upper 
panel) and 3a (lower panel) detections are displayed. The 
shadow plots show the stars with detected abundances, while the 
blanks also include stars with upper limits. 



common stars, our thorium abundance results are all in agree- 
ment with literatures values within the uncertainty as shown in 
Fig. |2j such consistence gives us the confidence that our results 
are reliable for studying the overall thorium abundances distri- 
bution of metal-poor stars. 

In Fig. [3] we plot thorium abundances and the ratios of 
thorium over europium against increasing metallicity, and we 
also plot results collected f rom previous stud i es in green, in- 
cluding field halo stars (e.g. .lCavrel et al.ll200lblHiU et al.ll2002t 
iFrebeletal .112 007: Roedere reTaLl l2009h . globular clusters M5 
and Ml 5 ('Vong et al. 2008a,l?), and a member of a nearby dwarf 
spheroidal galaxy (Aoki et al. 2007). Filled symbols are used to 
mark the common stars listed in |2] while due to different metal- 
licities adopted, they look not in pairs. In the lower panel a "zero- 



E 






] 






±- 





-1.5 

logE(Th)_THISWORK 



Fig. 2. Thorium abundances comparison for the common stars 
between this work and previous works, which are the averaged 
value if there are more than one results. Error bars are given and 
upper limits are marked with arrows. 



age" line caculated based on the initial ratio loffe(Th/Eu)„ = 
—0.33 given bv lSchatz et al.l (12002) is shown, and we will dis- 
cuss more about age estimations in section l4!4l Stars with both 
detected and upper limits of thorium abundances at the > Icr 
level are included. Europium and other element abundances in 
this work are from the analysis of Paper II. A summary of pre- 
vious results on thorium abundances is given in online Table |6] 

In the upper panel of Fig. [3] we see our sample provides de- 
tections which extend the range of derived thorium abundances 
to lower metallicities. At [Fe/H] < —1.8, a larger scatter can 
be seen, as expected according to the scenario of poor mixing 
in the very early universe. No clear trend with increasing metal- 
licity is seen, but the averaged detected loge(Th) is ~ 0.4 dex 
lower than those with metallicity above —1.8. Due to lack of 
derived Eu abundances, stars at the low metallicity end, e.g., 
HE 0353-6024 with the lowest metallicity of -3.17 in the 



4 



J. Ren et al.: The Hamburg/ESO R-process Enhanced Star survey (HERES) 



O r-I stars 

^ r-normal stars 



□ r— II stars 

O unclear type stars 



9Po 



oo 

o 



O 1 

oo,- 



c 



-2.5 
[Fe/H] 



o 



-2.5 
[Fe/H] 



Fig. 3. Plot of loge(Th) vs. [Fe/H] and log(Th/Eu) vs. [Fe/H]. 
Detected results from this work are in black, and the upper lim- 
its are in red. Results from other studies are in green. These 
filled symbols are the stars in common with previous works. 
See text for the definations of different subclasses. A dashed red 
line is plotted in the lower panel corresponding to the "zero-age" 
log(Th/Eu). Averaged error bars are given at the lower right cor- 
ner. 

Table 3. Definition of subclasses of metal-poor stars 



Class 


Constraints 


r-normal 


< [Eu/FeJ < 0.3 and [Ba/EuJ < 


r-I 


0.3 < [Eu/Fe] < 1.0 and [Ba/Eu] < 


r-II 


[Eu/Fc] > 1.0 and [Ba/Eu] < 


s 


[Ba/Fe] > 1.0 and [Ba/Eu] > 0.5 


S-II 


[Eu/Fc] > 1.0 and [Ba/Eu] > 0.5 


r/s 


0.0 < [Ba/Eu] < 0.5 



upper panel were lost in the lower panel. The distribution of 
log(Th/Eu) ratios spans from —0.86 to 0.50, while considering 
the uncertainty, most of them are consistently close to the "zero- 
age" line, supporting the universal r-process pattern in the early 
universe, except for only two stars: r-I star HE 0105—6141 and 
HE 1332—0309 with no clear type yet, which have much higher 
loge(Th/Eu) of 0.49 dex and 0.50 dex, respectively. They might 
belong to the so-called "actinide boost" star group, which we 
will discuss again in section l4!4l 



4.3. Correlation with otiier neutron capture elements 

Correlations between abundances of thorium and other neutron- 
capture elements are used to explore the nature of nucleosyn- 



thesis, and impose contraints on astrophysical and nuclear mod- 
els. Chemical abundance analysis suggest that for elements with 
Z = 52 ^ 76, the abundance pattern is consistent with the scaled 
solar system r-process abundance distribution. Recent studies 
have found relatively low abundance levels for light elements 
(Z < 50) compared with heavy elements, which suggest two 
distinct types of r-process events, a main r-process for the ele- 
ments above the second-peak elements, and a weak r-process for 
lighter neutron capture elements. The enhancement of actinides 
with respect to the rare earth elements were also noticed in some 
stars. 

In this work, we also explored the correlations between tho- 
rium and other neutron-capture elements. In Fig. 21 we plot the 
ratios of Th abundances against Ba, Sr and Y. There are larger 
scatters among the results of Th correlated with Ba and Sr, and 
usually the scatter comes from the unclear type stars, which 
indicates the existence of multiple nucleosynthesis prcesses. 
Scat ter of [Ba/Fe] and [Sr/Fe] were also re ported by other stud- 
ies jMcWilliaml 119981: iNorris et al.1 1200 Ih . All these elements 
follow similar abundance pattern as thorium as expected, and 
stars of the same enhancement type have quite consistent dis- 
tribution, supporting that they might form from the materials 
experienced similar chemical enrichments. Elemental ratios of 
log(Th/Ba), log(Th/Sr) and log(Th/Y) for almost all the stars are 
enhanced compared with the ratios in solar system. Two stars: 
HE 2219-0713 and HE 0353-6024 are particular interesting. 
The former star has more Th than Ba, and the latter one has 
more Th than Sr. As shown in Fig.|2l the Th II line is quite weak 
in both of the two stars, and even we can see some compara- 
ble "noise" or "features" unfitted, but since we failed to find the 
possible missing lines here, and our line list is believed to be 
well established, thus we suspect the S/N is overestimated for 
them, consequently leading to an overestimation of the detec- 
tion significance. Better quality data for the two stars are needed 
to confirm this results. 

In Fig. |5] the three chronometer ratios log(Th/Eu), 
log(Th/La) and log(Th/Nd) were plotted. Except the two stars 
with high enhancement of log(Th/Eu), all the other stars have 
consistent distribution with literature results, and the smaller 
scatter in the distributions, indicate the three age indicator will 
give consistent age estimations. 



4.4. Implications for age estimations and nuclear 
astrophysics 

The errors in Th abundance as well as the large scatter in the de- 
rived log(Th/Eu) ratios for metal-poor stars lead to significant 
uncertainties in age estimates. A 0.2 dex error in log(Th/Eu) 
will induce an error of 9.3 Gyr in age, which could be even 
more if considering the uncertainty in the predicted initial pro- 
duction ratios (PRs) (Frebel et al. 2007). In our case, the error 
of log(Th/Eu) is estimated to be > 0.25, which corresponds to 
an error of at least > 11.7 Gyr in age determinations. In addi- 
tion, it is known that, due to the so-called "actinide boost" phe- 
nomenon, the log(Th/Eu) chronometer pair fails to give mean- 
ingful estimates of age for some objects, e.g., CS 30306—132, 
CS 31078-018, CS 31082-001 and HE 1219-0312, all giv- 
ing negative ages. These objects exhibit higher abundances of 
thorium and uranium with respect to the lantha nides, which is 
currently not understood. iRoederer et alJ (12009 ^ examined the 
Pb and Th abundances in 27 r-process only stars, and suggest 
that, deviation from main r-process affect at most only the ele- 
ments beyong the third r-process peak elements, i.e., Pb, Th and 
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Fig. 4. loge(Th) plotted against neutron capture elements: Ba, Sr and Y. A diagonal 1:1 relation (black), and the solar ratio (red) 
are plotted for comparision. Error bars of the results from this work and the literature are given in the lower right of each plot. The 
symbols are as same as in Fig. [3] HE 2219—0713 has higher Th than Ba, and HE 0353— 6024 has higher Th than Sr. 
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Fig. 5. loge(Th) plotted a gainst neutron capture elements: Eu, La and Nd. A diagonal 1:1 relation (black), the initial ratio 
(green'l dSchatz et ani2002h and the expected ratio assuming an age of 13 Gyr (red) are plotted for comparision. Error bars are 
also given in the lower right of each plot. The symbols are as same as in Fig. [3] The two stars that have much higher Th than Eu are 
HE 0105-6141 and HE 1332-0309. 



U, while the reason for the low amount of Pb in CS 31082—001 
is not clear yet. Pb isotops and actinide elements measurements 
for more metal-poor stars are necessary for better understanding 
the nucleosynthesis in such actinide enhanced stars, and to ex- 
plore how common this "actinide boost" phenomenon may exist 
in metal-poor stars. 

Despite the large uncertainty in using log(Th/Eu) chronome- 
ter, it is meaningful to give the distribution of age estima- 
tions for so far the largest old population sample with avail- 
able log(Th/Eu). In Fig. |6] against [Fe/H], derived ages and 
the deviations of our expected initial production ratios (PRs), 
from the present theoretic initial PRs are plotted. We caculated 
the expected PRs assuming an age of 13 Gyr for all the stars, 
since these stars are expected to be very old based on their 
low metallicities. For age estimation we adopted log(PR) = 
—0.33 fro m Sclial/ e t al. (2002), which is in the middle of 
other PRs dSneden et a l. 2003; Cowan et all 120021: iFrebel et alj 
l2007h . Different PRs would only change the absolute age de- 



terminations by smalle amounts, i.e., '^2-4 Gyr dFrebel et al.l 
2007). Stars (green) from previous studies are also included 
for comparison. Ages from log(Th/Eu), together with that from 
lo g(Th/La ) and log(Th/Nd), if available, based on the initial PRs 
of lSchatz et al. (2002), are plotted in black, red and blue respec- 
tively in the middle panel. Lines corresponding to different ages 
or in the lower panel, the zero-deviation lines corresponding to 
different age chronometers are plotted. 

It can be seen from the plot that, the ages derived from 
log(Th/Eu) for those metal-poor stars in our sample are in gen- 
eral consistent with literature values within the uncertainty, ex- 
cept for HE 0105-6141 and HE 1332-0309, and there is no 
apparent trend with increasing metallicity. It's seen that the de- 
rived ages distribute around 13 Gyr and scatter below Gyr, 
with several "negative age" stars as previously reported from 
other works, and also in this work. In upper panel of Fig. [3] 
HE 0105-6141 and HE 1332-0309 are in the higher end of 
continuous loge(Th) distribution, but the ratio of log(Th/Eu) 
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Fig. 6. Upper Panel: plot of ages derived from log(Th/Eu) 
against [Fe/H]; Middle Panel: ages derived from different 
chronometers; Lower Panel: plot of the differences between our 
expected initial production ratios and the theoretic predictions 
log(Th/Eu)o against [Fe/H] assuming age of 13Gyr for all the 
stars. Symbols and colours are as same as in Fig. [3] and in the 
middle panel, we use black, red and blue symbols for results de- 
rived from log(Th/Eu), log(Th/La) and log(Th/Nd) respectively. 



reach as high as ^ 0.5 dex, ^ 0.5 dex higher than other stars. 
Thus, if considering the limitations of our data set, only these 
two stars are very likely to be "actinide boosted". 

In Fig. 121 the Thll line profile fits for stars above Icr 
detection are displayed. The observed spectra are in solid black 
lines, and the noise on the spectra were overplotted in dotted 
black lines. A range of 3 A fit region is highlighted in bright 
yellow. The red solid lines are the best fit synthetic profile and 
the red dotted lines are the profile fits with only blendings, from 
which it's easily seen how severe the Th II lines are blended. 
Star name, parameters, detection confidence m and to', and S/N 
are shown in the top of each plot. HE 1332—0309 has a low S/N 
of 31, and the observed line profile could't be reproduced well 
with present line list, which may due to low S/N or unresolved 
blendings in the Th II line region. HE 0105—6141 with higher 
S/N of 55, however also suffers from noise. We notice the 
blendings around Th II for the two stars are very light, which 
may also indicate that, some missing blending components 
might lead to an overestimation of Th abundance. Thus, better 
understanding of the blendings and higher quality data of 



these objects need to be obtained in the future to verify this 
preliminary result. 

For the common star CS 31082-001, we derivedlog(Th/Eu) = 
-0.39, log(Th/La) = -0.57, log(Th/Nd) = -1.02, which 
is con sistent with —0.22, —0.60, and —0.91 from ISchatz et al.l 
(l2002h within the un certainty. Althoug h we got a very consistent 
Th abundance with iHill et al.l (l2002h . we adopted a 0.15 dex 
higher Eu abundance from Paper II, which leads the offset in 
log(Th/Eu). We adopted the zer o-age log(Th/X)o (X is the stable 
r-process elements) ratios from lSchatz et al.l (l2002h . where they 
used solar abundances Xq instead of the r-process model 
predictions Xq, therefore, the discrepancies in the estimated 
ages from different log(Th/X) reflect deviations of the observed 
stellar abundances from a solar abundance pattern. As shown 
in the middle panel of Fig. |6l the deviations are tiny and 
similar for most stars, for which two or three chronometers 
are available, and log(Th/La) always gives the youngest ages, 
indicating compared with Eu and Nd, observed La abundances 
always have larger deviations from solar La abundance. For 
HE 0105-6141, only log(Th/Eu) and log(Th/Nd) are available, 
the latter one gives a much older age, which may suggest the 
high log(Th/Eu) in this stars is due to an abnormal low Eu 
abundance, while so far no other studies give Eu abundance for 
this star. 

From the bottom panel of Fig. |6] it shows if assuming a 
consistent age of 13 Gyr for these metal-poor stars, their intial 
production ratios for most of the stars seems needed to be 
increased by a factor up to 4, and even more if the two "actinide 
boost" stars are confirmed. 

In order to explore more about these two unique stars, in Fig. 
[8] we plot their abundance patterns, compared with solar r- and 
s-process patterns scaled to Eu and Ba (or Sr) respectively. For 
HE 1332—0309, because no Ba abundance has been obtained, 
we scaled the solar s-process pattern to Sr. As comparison, a 
thorium normal r-II star CS22892— 052, a previous reported "ac- 
tinide boost" r-II star CS 31082—001 are also plotted. It can be 
seen from the plot, except for Th, the available neutron-capture 
abundances of CS22892— 052 follow the Solar pattern quite 
well, and the same is for CS 31082—001, but for both stars, we 
underestimated their log(Th/Eu) ratios. The situation is different 
for HE 0105—6141, for which the Th abundance was highly 
enhanced, but the other available neutron-capture elements 
follow the Solar r-process pattern. For HE 1332—0309, Barium 
abundance is not available. Beside Th, Nd, Sr, Zr, and Y are also 
enhanced to some degree compared with solar r-process pattern. 



5. Conclusion 

Using the "snapshot" spectra from HERES survey, we derived 
the thorium abundances of 77 metal-poor stars, 17 of which have 
detected Th abundance, while for the rest, only upper limits are 
available. Thorium abundances cover a wide range of about 4.0 
dex, and a scatter exists in the distribution of log(Th/Eu) ratios, 
supporting previous studies suggesting that the r-process is not 
universal. For the five common stars, our results are in good 
agreement with previous studies, which gives us the confidence 
to present a reliable Th distribution and to discuss r-process pat- 
tern for such a large sample of metal-poor stars. 

With avaliabe abundances, we explored the correlation be- 
tween Th and other r-process elements. We confirmed the rel- 
atively large scatter in log(Th/Ba) and log(Th/Sr) distributions, 
and found that a better consistence exists within r-process stars. 
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which confirms that these stars formed from the gas experienced 
similar nucleosynthesis. 

Using log(Th/Eu) as the chronometer, we derived the ages 
for r-process metal-poor stars, for which both Th and Eu are 
available. Two stars might be "actinide boosted", considering 
an error to 11.7 Gyr in our age estimation. These stars might 
have experienced very different chemical enrichments during 
their formation and evolution. At present the only possible ex- 
planation would be that, the r-process elements in these stars 
were implanted long after the formation of the star, still hold- 
ing the principle of a universal r-process pattern. It suggests that 
"actinide boost" might not be a rare phenomenon, thus question- 
ing the reliability in using log(Th/Xstabic) for cosmochronom- 
etry to derive stellar ages. For this kind of stars, U abundances 
are needed to confirm the "actinide boost" feature and log(Th/U) 
can be used to give more accurate age estimations. 

However, it should be cautioned again that, la detection sig- 
nificance was adopted in this work, due to the weakness of Th 
II line and the limited quality of our present data, which may 
lead to significant uncertainty in the thorium abundance. Thus, 
the detection of Th in stars without enhancement of other n- 
capture elements and the "actinide boost" phenomenon in some 
stars are not confirmed conclusion. Better quality spectra are 
needed for higher precision Th and other r-process heavy ele- 
ments abundance determinations. Thus here we would not like 
to give any conclusion, but rather to draw attention on those 



very interesting objects from this work, i.e., HE 2219—0713, 
HE 0305-6024, HE 0105-6141 and HE 1332-0309. Actually 
we ourselves have already submitted a proposal for the obser- 
vation time using ESO VLT/UVES for obtaining better quality 
data of these objects, if some of them can be confirmed, it will 
be a very important discovery. 
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Table 5. Summary of abundance results for stars with Th abundances from this work. Except for Th, other element abundances and the stellar 
parameters are from Paper II. The errors are discussed in the text, m and m' are the detection levels for Thll 4019.129 line profile fit with and 
without blendings. 
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Table 5. Continued. 



Name 


7c ff 


log5 


[Fe/H] 




loge(Th) 


[Th/Fe] 


[Eu/Fe] [BaTEu] log(Th/Eu) [C/Fe] 


m(CT) 


m'ia) 




[K] 


[cm s~^] 




[kms-i] 




















HE 1431-2142 


6137 


4.10 


-2.60 


1.55 


< 


-0.20 


< 2.31 


- 


- 


- 


0.48 


0.60 


0.23 


HE 2133-1432 


5716 


3.46 


-2.02 


1.46 


< 


-0.82 


< 1.11 


- 


- 


- 


0.12 


1.10 


0.60 


HE 2134+0001 


5257 


3.00 


-2.22 


1.55 


< 


-1.36 


< 0.77 


0.47 


-0.89 


< -0.12 


0.20 


1.17 


0.52 


HE 2151-2858 


5598 


4.14 


-2.38 


1.39 


< 


-0.43 


< 1.86 


- 


- 


- 


0.10 


0.68 


0.51 


HE 2153-2719 


4898 


2.01 


-2.49 


1.80 


< 


-2.81 


< -0.41 


0.21 


-0.90 


< -1.04 


0.12 


0.95 


0.24 


HE2158-3112 


4843 


1.85 


-2.75 


1.85 


< 


-2.81 


< -0.15 


0.02 


-0.86 


< 0.60 


-0.04 


1.31 


0.42 


HE 2219-0713 


4789 


1.68 


-2.91 


1.64 




-1.97 


0.85 


- 


- 


- 


-0.17 


1.35 


0.97 


HE 2224+0143 


5198 


2.66 


-2.58 


1.67 




-1.47 


1.02 


1.05 


-0.46 


-0.45 


0.35 


1.72 


0.99 


HE 2224-4103 


5074 


2.32 


-2.64 


1.75 


< 


-1.71 


< 0.84 


- 


- 


- 


0.23 


0.90 


0.45 


HE 2229-4153 


5138 


2.47 


-2.62 


1.79 


< 


-1.91 


< 0.62 


0.45 


-0.73 


< -0.28 


0.37 


0.79 


0.20 


HE 2234-0521 


5332 


3.15 


-2.78 


1.39 


< 


-0.93 


< 1.76 


- 


- 


- 


0.36 


0.27 


0.07 


HE 2247-3705 


5366 


3.04 


-2.27 


1.48 


< 


-1.39 


< 0.79 


- 


- 


- 


0.36 


0.76 


0.27 


HE 2250-2132 


5705 


3.69 


-2.22 


1.43 


< 


-0.74 


< 1.39 








0.41 


0.89 


0.43 


HE 2252-4157 


5090 


2.87 


-1.93 


1.55 




-1.32 


0.52 


0.53 


-0.24 


-0.43 


-0.15 


2.98 


1.34 


HE 2259-3407 


6266 


4.32 


-2.29 


1.37 


< 0.00 


< 2.2 








0.41 


0.39 


0.16 


HE 2311+0129 


5188 


2.65 


-2.78 


1.54 




-1.32 


1.37 








0.33 


1.24 


1.05 


HE 2326+0038 


5145 


2.51 


-2.77 


1.62 


< 


-1.59 


< 1.09 








0.23 


0.51 


0.11 


HE 2327-5642 


5048 


2.22 


-2.95 


1.69 




-1.45 


1.41 


1.22 


-0.56 


-0.25 


0.43 


1.79 


1.35 


HE 2338-1618 


5515 


3.38 


-2.65 


1.43 


< 


-0.76 


< 1.80 








0.47 


0.88 


0.56 


HE 2345-1919 


5617 


4.46 


-2.46 


1.47 


< 


-0.47 


1.90 








0.24 


0.97 


0.60 
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Table 6. Thorium and europium abundances from previous works (From the summary of SAGA database and lRoederer et alJ2009l) . 
Both detected values and upper limits are included. 



name 


[Fe/H] 




loge(Th) 




loge(Eu) 


log(Th/Eu) 


"Age" 
[Gyr] 


Ref. 


BD+08.2856 


-2 


10 


-1 


78 







10 


-1 


16 







06 




-0 


62 


13.5 


1,2 


BD+17.3248 


-2 


23 


-1 


18 







10 


-0 


67 







05 




-0 


51 


8.4 


3 


BD+04.2621 


-2 


50 




< 




3 


12 


-2 


61 







23 


< 


-0 


51 


> 8.4 


1,2 


BD- 18.5550 


-3 


03 




< 




3 


10 


-2 


79 


± 





22 


< 


-0 


31 


> -0.9 


1,2 


CS22892-052 


-2 


92 


-1 


42 


± 





15 


-0 


86 







02 




-0 


56 


10.7 


7 


CS29491-069 


-2 


60 


-1 


43 


± 





22 


-1 


03 


± 





10 




-0 


40 


3.3 


3 


CS29497-004 


-2 


81 


-0 


96 







15 


-0 


45 







20 




-0 


28 


-2.3 


4,5 


CS30306-132 


-2 


40 


-1 


12 


± 





15 


-1 


14 


± 





25 







02 


-16.3 


7 


CS31078-018 


-2 


84 


-1 


35 


± 





25 


-1 


17 


± 





17 




-0 


18 


-7.0 


8 


CS31082-001 


-2 


90 


-0 


98 


± 





05 


-0 


76 


± 





11 




-0 


22 


-5.1 


8,9,10 


HD108317 


-2 


18 


-1 


84 







20 


-1 


32 


± 





05 




-0 


52 


8.9 


3 


HD108577 


-2 


36 


-2 


11 


± 





14 


-1 


48 


± 





02 




-0 


63 


14.0 


1,2 


HD110184 


-2 


52 


-2 


50 







15 


-1 


91 







05 




-0 


59 


12.1 


7 


HD 11 5444 


-2 


85 


-1 


97 


± 





15 


-1 


64 


± 





03 




-0 


33 


0.0 


3 


HD122563 


-2 


72 




< 




2 


43 


-2 


75 







11 




< 


32 


> -30.3 


3 


HD122956 


-1 


95 


-1 


50 







17 


-0 


94 







07 




-0 


56 


10.7 


3 


HD 1265 87 


-2 


93 




< 




2 


39 


-1 


97 







06 


< 


-0 


42 


> 4.2 


3 


HD128279 


-2 


00 




< 




2 


00 


-1 


57 







06 


< 


-0 


43 


> 4.7 


2,18 


HD175305 


-1 


48 


-0 


76 


± 





15 


-0 


36 







07 




-0 


40 


3.3 


3 


HD 186478 


-2 


50 


-1 


85 


± 





15 


-1 


34 


± 





06 




-0 


51 


8.4 


7 


HD204543 


-1 


87 


-1 


68 


± 





14 


-1 


05 


± 





07 




-0 


63 


14.0 


6 


HD221170 


-2 


20 


-1 


46 







05 


-0 


86 







07 




— 


60 


12.6 


11 


HD6268 


-2 


40 


-1 


93 


± 





10 


-1 


56 


± 





03 




-0 


37 


1.9 


3 


HD74462 


-1 


52 


-0 


94 


± 





13 


-0 


50 


± 





09 




-0 


44 


5.1 


3 


HD88609 


-3 


07 




< 




2 


65 


-2 


89 


± 





12 




< 


24 


> -27 


19 


HE0338-3945 


—2 


42 











23 





02 







17 




< 


21 


> -25 


6 


HE1219-0312 


-2 


97 


-1 


29 


± 





14 


-1 


06 


± 





10 




-0 


23 


-4.7 


3 


HE1523-0901 


-2 


95 


-1 


20 


± 





05 


-0 


62 


± 





05 




-0 


58 


11.7 


12 


HE2 148 -1247 


-2 


50 




< 







50 





17 


± 





10 


< 


-0 


67 


> 15.9 


20 


HE2327-5642 


-2 


78 


-1 


67 


± 





21 


-1 


29 


± 





07 




-0 


38 


2.3 


13 


M15 K341 


-2 


32 


-1 


51 


± 





10 


-0 


88 


± 





09 




-0 


63 


14.0 


14 


M15 K462 


-2 


25 


-1 


30 


± 





10 


-0 


61 


± 





09 




-0 


69 


16.8 


14 


M15 K583 


-2 


34 


-1 


70 


± 





10 


-1 


24 


± 





09 




-0 


46 


6.1 


14 


M5 IV-81 
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